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We develop a general theory of fermion liquids in spatial dimensions greater than one. The prin- 
cipal method, bosonization, is applied to the cases of short and long range longitudinal interactions, 
and to transverse gauge interactions. All the correlation functions of the system may be obtained 
with the use of a generating functional. Short-range and Coulomb interactions do not destroy the 
Landau Fermi fixed point. Novel fixed points are found, however, in the cases of a super-long 
range longitudinal interaction in two dimensions and transverse gauge interactions in two and three 
spatial dimensions. We consider in some detail the 2-|-l-dimensional problem of a Chern-Simons 
gauge action combined with a longitudinal two-body interaction ^(q) oc |q|^~^ which controls the 
density, and hence gauge, fiuctuations. For y < we find that the gauge interaction is irrelevant 
and the Landau fixed point is stable, while for j/ > the interaction is relevant and the fixed point 
cannot be accessed by bosonization. Of special importance is the case j/ = (Coulomb interaction) 
which describes the Halperin-Lee-Read theory of the half-filled Landau level. We obtain the full 
quasiparticle propagator which is of a marginal Fermi liquid form. Using Ward Identities, we show 
that neither the inclusion of nonlinear terms in the fermion dispersion, nor vertex corrections, alters 
our results: the fixed point is accessible by bosonization. As the two-point fermion Green's function 
is not gauge invariant, we also investigate the gauge-invariant density response function. Near mo- 
mentum Q = 2fc_F, in addition to the Kohn anomaly we find singular behavior. In Appendices we 
present a numerical calculation of the spectral function for a Fermi liquid with Landau parameter 
/o 7^ 0. We also show how Kohn's theorem is satisfied within the bosonization framework. 
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I. INTRODUCTION 

The peculiar nature of electrons in a half-filled Landau level, and anomalies in the normal state transport properties 
of the high-temperature cuprate superconductors, have stimulated great interest in systems of strongly interactiDg 
fermions moving in spatial diiaensions greater than one. In addition to methods such as mean-field theorylil, RPAld'c 
and 1/N perturbation theogya-tl, a number of non-perturbative theoretical tools such as the renormalizatioa group,QE 
the Eikonal approximation^ and Ward identitieaj have been used to study fermion liquids. BosonizationCTllJ is another 
powerful tool for understanding the effects of both regular and singular interactions on fermion liquids as the realm 
in which it is applicable, for low-energy excitations near the Fermi surface, is precisely where these interactions are 
most important. Bosonization does not rely upon a Fermi liquid form for the quasiparticle propagator; consequently 
non-trivial zero temperature quantum critical fixed points are accessible. In an appropriate limit, bosonization reduces 
complicated many-body Hamiltonians to quadratic forms which possess infinite U{1) symmetry and which can be 
diagonalized exactly. 

A general theory of fermion liquids in two and three spatial dimensions requires the consideration of a variety of 
Fermi surface geometries and interactions. In this paper we study only the simplest Fermi surface geometries: a circle 
in spatial dimension Z? = 2 or a sphere in D = 3. We assume that the symmetries of these surfaces do not break 
spontaneously when the interactions between fermions are turned on. The bosonization formalism applies, without 
significant modification, to more complicated Fermi surfaces provided there are no Van Hove singularities and nesting 
does not occur. However, the added technical complications of describing these more complicated surfaces would make 
detailed discussion of their properties cumbersome and, given this caveat, would not lead to any qualitative changes 
in the physics. We do consider a variety of two-body interactions between the fermions. In accord with expectations, 
we find, barring superconducting instabilities, that short-range and Coulomb interactions do not destroy the Landau 
fixed point. Novel fixed points are found, however, in the cases of super-long range interactions and transverse gauge 
interactions. The latter case is particularly significant as gauge fields play an. important role in the physics of the 
half-filled Landau levelEl and may be important in the cuprate superconductorsESEj. The effects of transverse gauge 
fieldS|r-svhich are not screened, are particularly interesting as perturbative expansions in the coupling constant break 
downEJ. 
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Earlier work on multidimensional bosonization showed that it is an accurate calculational framework for understand- 
ing Landau Fermi liquidsLrEZI. Here we address the question of whether the approximations inherent in bosonization 
are adequate for the method to apply to a wider range of fermion liquids with singular or transverse gauge inter- 
actions. Three questions in particular must be answered in order to check that bosonization does not break down. 
First, the assumed current algebra must be consistent: it is derived for the free system, and it must continue to hold 
when interactions are turned on. Second, we examine whether or not the neglect of quadratic terms in the fermion 
quasiparticle spectrum due to Fermi surface curvature is an innocuous approximation. Finally, we need to check 
whether certain vertex corrections are significant by returning to the fermion basis to examine their contribution. 
In this paper, we examine systems which meet all three of the above conditions. The fixed points of these systems 
therefore can be accessed and studied via bosonization. 

The outline of this paper is as follows. In Sec. ^ we develop a bosonization formalism applicable in any spatial 
dimension which describes both longitudinal and transverse gauge interactions. A generating functional permits the 
easy computation of 2N-point boson and fermion correlation functions. Interactions of interest are simply special cases 
of this general formalism. In Sec. Ill, we demonstrate how to incorporate any two-body density-density interaction 
by introducing longitudinal gauge fields which mediate the interaction. As examples, we reproduce the Luttinger 
liquid in one spatial dimension and solve the problem of a single, but general. Landau parameter /„ in two and three 
spatial dimensions. In Appendix A we use the numerical method of fast-Fourier-transforms (FFT) to extract the 
spectral function of a two-dimensional Fermi liquid with Landau parameter /o. And in Appendix B we present details 
of our solution for more complicated Landau parameters. Then we turn to the problem of long-range longitudinal 
interactions. We show that screening of the Coulomb interaction occurs naturally within the bosonization framework. 
By considering fermions with spin, we show that spin-charge separation does not occur in D > 1. In the case of 
super-long range interactions in D = 2, however, the Landau fixed point is destroyed by the emission of plasmons by 
the quasiparticles. We show that bosonization accurately describes this breakdown. 

The fact that bosonization can access non-Fermi liquid fixed points, both in one and two spatial dimensions, 
motivates Sec. [V, in which we study a more physically relevant example of the breakdown o£ Fermi liquid behavior: 
the Chern-fSimons theory of half-filled Landau level put forward by Halperin, Lee, and ReadlJ (HLR) and Kalmeyer 
and Zhangp2l. To investigate the range of validity of bosonization, we examine the general case of an interaction of the 
form V{q) = |q|^^^ with — 1 < y < 1 when the inverse gauge propagator is given by i7|ijj|/|q| — xlll^^^- Aided by the 
Ward Identity approach, we find that for y < the gauge interaction is irrelevant and the Landau fixed point is stable, 
while for y > the interaction is relevant and the fixed point cannot be directly studied by bosonization. We discuss 
in detail the special case y = corresponding to the Coulomb interaction. In this case a marginaLEermi liquid (MFL) 
fixed point controls the low energy physics. (An brief account of this work has appeared in printllj.) However, as the 
two-point Green's function is not gauge invariant, it is not directly related to experimentally observable quantities. 
The density response function, on the other hand, is gauge invariant and in Sec. \4we study its behavior for Q « 
and Q Ki 2kF- In the former case our results are identical to those found inJlPAH while in the latter case we find 
non-analytic structure stronger than that found by Altshuler, loffe and MillisH. In Sec. VI, we apply the formalism 
developed in Sees. || and IV to the case of the physical Maxwell electromagnetic interaction in three dimensions 
to obtain the quasiparticle Green's function in the Coulomb gauge. Although the Landau fixed point is unstable at 
extremely low energies, in accord with expectations it contro ls th e behavior at physically relevant temperatures. A 
brief summary of our work and conclusions are found, in Sec. VIL In Appendix C, we include the Landau parameter 
/i in the HLR theory to show that Kohn's theoremcj is recovered automatically: the bare cyclotron frequency, not a 
renormalized value, appears in the collective mode spectrum. 



II. BOSONIZATION 



In this section we formulate the problem of an interacting fermion liquid, making use of a bosonization method 
applicable in any spatial dimension to study both longitudinal and transverse gauge interactions. A generating 
functional permits the ready computation of 2N-point boson and fermion correlation functions. 

We start from the bare Hamiltonian for fermions interacting with non-compact longitudinal and transverse gauge 
fields A'^ in _D-spatial dimensions. For simplicity, we consider spinless fermions and spherical Fermi surfaces only. In 
the Coulomb gauge, V • A = 0, and when % is set equal to one, 

if = y'd^xct(x)[t^^--^-Ao]c(x) + i l^^xd^y F(x-y) ct(x) ct(y) c(y) c(x) (1) 

where ek = k^/2m, and l^(x — y) represents Coulomb or other two-body inteeactions. Next we make use of the 
renormalization group to integrate out the high-energy Fermi degrees of freedomQ. The resulting low-energy effective 
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theory is expressed in terms of quasiparticles ip^ which obey canonical anticommutation relations and which are 
related to the bare fermion operators by 



= Ck (2) 

for momenta k which are restricted to a narrow shell of thickness A around the Fermi surface: fc^— A/2 < |k| < kF+X/2. 
The wavefunction renormalization factor rescales the discontinuity in the quasiparticle occupancy at the Fermi 
surface back to one. To be explicit, consider the partition function 

Z = J VA^" VipV ^ e' ^' (3) 

where we have separated the fermion fields into low and high energy components, tpi^) ^J^d //'(x) respectively. 
Integrating out the high energy /p fields, which contain all of the degrees of freedom except for the narrow shell 
around the Fermi surface, we obtain the effective action: 

S[A,tp]^ fd^xdt [r{^)(idt + AoM^) + ^ t AftPi^)] 

,1 2m* 

(4) 



We use the notation q= (oj, q). Renormalization of the Fermi velocity vp ^ v*p is expressed through the appearance 
of the effective mass m* . The last term in Eq. (^, which is generated by integrating out the fields deep inside 
the Fermi sea contributes to the total diamagnetic term, — (py/2m*)A(g)-A(— g); pf is the fermion number density. 
Irrelevant operators consisting of higher powers and/or derivatives of the gauge and Fermi degrees of freedom are also 
generated. These |«Berators may be safely neglected in the low-energy limit assuming that the anomalous dimension 
of the gauge fieldsOla is not too different from their engineering dimension. To check that this-is the case, the effective 
action of the gauge fields must be obtained by integrating out the fermion degrees of freedomQO, and examined for the 
case of interest. The effects of two-body fermion interactions which in general include both the long range Coulomb 
interaction and short range Fermi liquid type interactions will be considered later in thcr-discussion. 
To bosonize the effective action Eq.(^) we introduce the coarse grained charge currentll3 

J(S; q) = E ^(S; k + q) ^(S; k) {^1+^ Vk - <^^,o "k} • (5) 

k 

Here S labels a patch [S = (0, (f) in three dimensions] on the Fermi surface with momentum kg, and 0(S; k) = 1 if k 
lies inside a squat box centered on S with height A in the radial (energy) direction and area A^~^ along the Fermi 
surface, and equals zero otherwise. These two scales must be small in the sense fcj? ;» A ^ A > |q|: we satisfy these 
limits by setting A = kp/N and A = kp/N" where < a < 1 and N oo. The currents so defined satisfy the U{1) 
current algebra 

[J(S; q), J(T; p)] = <5^t'Cp,o" ^ ' ■ (6) 

Here f2 = A^~^{L/2Tr)^ is the number of states in the squat box divided by A and fis is a unit normal to the Fermi 
surface at S. The charge currents also have a bosonic representation; if we set 

J(S;x) = Vi^ns • V(/)(S;x) (7) 

where the bosonic charge fields satisfy the canonical commutation relations. 



_( 

= 0; \x^-y^\A:^l (8) 



where e{x) = 1 for x > and equals —1 otherwise. Eq. (g) follows immediately from the commutation relations. The 
key formula of the bosonization procedure expresses the fermion quasiparticle field ip in terms of the boson fields as 



^(S; X, i) = ^7 V -^""'^ exp{z^</.(S; x, t)} 0(8) . (9) 
i/y V a It 
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Here V is the volume of the system and a = 1/A is an uhraviolet cutoff. 0{S) is an ordering operator introducedl3'E^I 
to maintain Fermi statistics in the angular direction along the Fermi surface. Anticommuting statistics are obeyed 
automatically in the direction normal to the Fermi surface, just as in one-dimensional bosonization. It is now an 
interesting exercise to use this representation of the fermion field operators to verify that 

J(S; x) = y lim {V'^(S; x + eng) ^{S; x) - (V/(S; x + eng) ^{S; x)) } 

= ^/4^ns • V(/)(S;x) . (10) 

The currents are invariant under U{1) phase rotations through an angle /3(S) in each patch: iIj{S; x) e^^^^^ V'lS; x). 
This infinite C/(l) symmetry reflects the fact that the current operator in a given patch does not scatter quasiparticles 
outside of that patch. The effective action Eq. (^ now can be expressed in terms of the charge currents J and 
consequently reflects the underlying infinite U{1) symmetry of the fixed point: 

S[A,a]^J2 E / J (^-fFq-fts)a*(S;g)a(S;g) 

+ / |^{E'^(S;(z)[Ao(-g)+v^A(-g)]--^A(g).A(-g)} . (11) 

Here Vg = ks/m* is the renormalized Fermi velocity vector at patch S and the charge currents are related to the 
canonical boson operators a and by: 

J(S; q) = |ns-q| [a(S; q) e{hs ■ q) + a^(S; -q) 0(-ns • q)] , (12) 

where 

[a(S;q), a+(T;p)] =5^^' C (13) 

and d{x) = 1 if a; > and is zero otherwise. Exchange scattering, which involves high momentum gauge fields, has 
not been included in the action, as it is not singular in the low-momentum limit. The action Eq. ( |ll| ) must be 
supplemented by the bare gauge action which we take to have the quadratic form: 

- ^ / (0^ / S ^'^'^^ ^^^"'^^ ■ 

We fix the gauge with the supplemental condition V • A = (ic. the Coulomb gauge). Particular examples considered 
in this paper are the physical electromagnetic gauge action in = 3, and in D = 2 the Chern-Simons action, which 
breaks time-reversal and parity, but not charge-conjugation, symmetries. 

We proceed with the evaluation of the boson correlation function which can be carried out exactly. Rather than 
integrate out the gauge fields to obtain an effective current-current interaction, we introduce a generating functional for 
the boson correlation function and first integrate out the boson fields a and a* . In this way we avoid the intermediate 
step of summing an infinite perturbation series in the effective interaction to obtain the boson propagatorllil. We 
construct the generating functional by coupling fields f and ^* to the boson fields a* and a. 

Z[£.,C] ^ J T^'A" Va Va* exp{i{S[A,a] + S^[A])} 

xexp{*E E 1^ K(S;g)a*(S;g)+r(S;<z)a(S;q)]} . (15) 

S q.q ns>0 

Here the prime over the gauge measure indicates that the integration respects the Coulomb gauge fixing condition. 
Then by completing the square we integrate out the boson fields in each patch S. It is important to note that the 
entire Fermi surface participates: patches at every point on the Fermi surface contribute to the effective action. 

d^q f duo 6{(\ ■ ns) 



Z[LC]^^ J V'A^ exp{z 5G[A]}exp| J ^ J 



{2tt)^ J 271 lu- v*pCi ■ ns 
X ^flq-hs K(S;g) (Ao(-g) + v^A(-q))+r(S;g) (^0(9) + A(g))] 

-i-ye(S;g) r(S;g)ll , (16) 
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where S'g[^] is an effective gauge action, wliicli since tfie bosons were coupled linearly to the gauge fields is a quadratic 
form, 



The inverse of the gauge propagator satisfies the equation 

K^,{q) Af^iq) A'^i^q) = K^q) A^{q) A'' {-q) + (q) Ao{q) Ao{-q)+X^{q) A{q) ■ A{~q) 



(17) 



(18) 



where x° and are the longitudinal and transverse susceptibilities which for the cases of perfectly circular or 
spherical Fermi surfaces, are given by 



x"{q)^N*mi 



1) \x\ 



a 



£> = 3 , 



D 



(19) 



and 



X^{q) = v*/ N*{0)[- x^ + e{x^ -1) \x\Vx^ -l + ieil-x"^) \x\VT^] 
= lv*^'N*{0)[-x'-^{l-x')ln{^^)] ; D = 3, 



D 



(20) 



where x = uj/{vp |q|), and A^*(0), which is the quasiparticle density of states at the Fermi surface, equals m*/27r for 
spinless fermions in D = 2, and m*fci?/27r^ in _D = 3. The integral over the gauge fields can be carried out now 
giving the generating functional as an explicit function of ^ and ^* . The boson correlation functions are obtained by 
differentiating the logarithm of the generating functional with respect to ^ and In particular the boson propagator 
defined by 



{a{S;q) a\S;q)) 



SCiS;q) (5e(S;g) 
is given in terms of the gauge propagator D^y{q) = [if (g)^^] ^ as 



{a{S-q) at(S;<7)) 



uj - v*pC\ • ns + 177 sgn(a;) {2tt) 



. AP-^ ^ i^^,(q,c^) e^(S;g) e'^(S;-g) 
I , ^ q ■ ns - 



v*p<i -hs+irj sgn(ti;)]2 



(21) 



(22) 



where e^{S\q) is the D+l-vector (l,Vg). Note that the velocity of the interacting bosons differs from that of free 
bosons by at most of order v*p {A/kp)-^^^ . Therefore velocity renormalization of the bosons due to interactions among 
the final shell of excitations with energy less than VpX is insignificant in the A ^ limit except in the special case 
of one spatial dimension. We will return to this observation below when we address the question of whether or not 
spin-charge separation occurs in two or higher spatial dimensions. 

Now the bosonization formula Eq. (^ can be used to express the space-time fermion Green's function in terms of 
the Fourier transform of the boson correlation function. In general terms the Green's function for fermions in a single 
patch on the Fermi surface is given by 



Gj^(S;x,<) = Ae^ksxg^p 
V a 



d^q 



2n ^ 



{2tt 



ns • q 



{a{S-q) a\S-q)) 



(23) 



which becomes, upon using Eq. (|22| ) and dropping the prefactor of Vl/V for brevity. 



Gi.(S;x,t) 



x • fis - v*pt 



oiks-x 



x • ns 



V*pt 



exp 



exp 



{2ti)D 
<5Gs(S;x, t) 



2^ 



gi(q-x-(^t) 



1] 



D^M,^) e''(S;g) e^(S;-g) 



(24) 



for x^ = x X fis = 0. The condition x_l = can be replaced by |x_l|A < 1 provided we control the logarithmic infrared 
divergence of the free boson correlation function by placing the system inside a large box. The weak divergence, which 
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is a consequence of treating the Fermi surface as locally flat inside each patch, presents no real difficulties in practice, 
and will be ignored in the following analysis. For |xx|A ^ 1 the fermion Green's function vanishes since bosons 
separated by large xj^ are uncorrelated. In the second line of Eq. (|2j) for convenience we have introduced the 
notation SGb for the modified Fourier transform of the additive correction to the free boson propagator due to 



interactions. 



d^q fduj 

A/2 



5Gb{S; ^,t)=t I J — [e'(i— *) - 1] 6Gb{S; q, cj) , (25) 



where we can identify SGb{S; q, to) as: 



dGij(S;q,w) = ; — — tTiT • (26) 

In Eq. (I25I), the momentum integral over q ranges over (—A/2, A/2) in each direction. In contrast, in the free part 
of the boson and fermion propagators, momenta in the perpendicular directions range over the much larger interval 
(—A/2, A/2). The reason for the difference lies in the fact that interactions which couple different patches on the 
Fermi surface must be cutoff at momentum A <SC A because in general the Fermi surface normal vector points in a 
different direction for each patch. Therefore, as the patches are squat pillboxes with dimensions A-^^^ x A, only modes 
with |q| < A are permitted by the geometry of the construction to couple patches together. 

To complete the program of computing the fermion two point function, the Green's function in Eq. ( p^ is summed 
over all patches on the Fermi surface: 

Gf(x,i) Gf(S;x,t) . (27) 

s 

The prime indicates that the sum is only over patches for which |x x-nsj A < 1. The expression for the fermion Green's 
function given here coincides with that found by Castellani et al.u who make use of Ward Identities derived in the 
fermion basis to find Eq. (E4). Why these two different approaches should give the same result and the consequent 



implications are discussed in detail in Sec. IV D 



III. LONGITUDINAL INTERACTIONS 



In this section we show how the generating functional may be used to determine the correlation functions of a 
fermion liquid interacting via two-body forces in two important special cases, for short-range interactions of Landau 
Fermi liquid type, and for Coulomb or longer range longitudinal interactions. The procedure parallels that given in 
the previous section as we introduce fictitious gauge fields, which are linearly coupled to the bosons, to mediate the 
interaction. First, to illustrate the non-perturbative nature of the method, we determine the propagator of a fermion 
liquid in one spatial dimension, which has a Luttinger liquid form. In higher dimensions, we extend the discussion to 
general short-range interactions which as expected do not destroy the Fermi liquid state. The quasiparticle propagator 
for a system of fermions interacting via a short-range interaction described by a single Landau parameter /o had been 
derived in a previous worktiJ by a different technique. Here we give a detailed derivation of the propagator for the 
nontrivial case of a fermion liquid with an /i interaction. This result will be used later in the derivation of Kohn's 
theorem which is applied in the discussion of the half-filled Landau level studied in the next section. The derivation 
for a general interaction specified by more than one Landau parameter and a proof of Kohn's theorem are given in 
Appendices B and C respectively. Finally, we show that, due to screening, the ground state of a fermion system 
interacting via the Coulomb interaction is a Landau Fermi liquid in two and three dimensions. For a super-long-range 
interaction in D = 2, however, the emission of plasmons by quasiparticles destroys the Fermi liquid fixed point. 

In one dimension, even a short-range interaction drives the fermion system to a different fixed point, which is a 
Luttinger liquid. The formalism developed in Sec. ^ can be used to diagonalize the interacting fermion problem 
exactly and to access the Luttinger liquid. As an example we consider an interaction which makes a contribution to 
the action of the form 

^-=-y/SI ^ JiS;<l) J(T;-q) (28) 

where L and R specify the left and right Fermi points. We introduce a fictitious field, A^^q), which generates this 
interaction when integrated out of the partition function. 
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s' 



2-K J 27r [ 2/( 



■Ao[q) Ao(-g) + [J(i; q) + J{R; q)] Ao(-g) 



(29) 



Now the total action has the form discussed in Sec. || and we follow the procedure given there to determine the boson 
correlation function. We find the propagator for the right hand patch to be given by 



{a{R;q) a^{R;q)) = 



^l-Fi/A^uj + qv'p^l-F^/A- 



(30) 

where v'p — vp + /o/(27r) is a renormalized Fermi velocity due to intrapatch interactions and Fg — fo/iTTv'p). If we 
now define r] such that 



cosh'^ V = -^[^ 



the boson propagator can be written as 



{a{R; q) a^R; g)) = ^| cosh^ j] 



u;~qv'p^l-F-^/4 



sinh^ 7] 



LO + qv'p^l-F^/A- 



(31) 



(32) 



where tanh277 = Fq. The fermion propagator found from Eq. ( P3| ) is non-analytic in the coupling constant /o and 
cannot be obtained by perturbing in the coupling constant, 



Gp{R; X, t) oc 



X — V*pt 



(33) 



where the anomalous exponent a = sinh^ rj and the Fermi velocity Vp = v'p yj\ — F^ /A. The repilt is exact and 
identical to that found by diagonalizing the bosonic Hamiltonian with a Bogoliubov transformationlHJ. 

A&_a second example, we consider a fermion liquid in two dimensions with short-range interactions. In earlier 
workO we discussed, by a different but equivalent approach, the effect of an interaction which could be parameterized 
by a single Landau parameter, /q. It was shown that a perturbative expansion in /o gave for the imaginary part of 
the self-energy a term proportional to f^uP' \n\uj\ in two dimensions and /qW^ in three dimensions as in Fermi liquid 
theory. A non-perturbative treatment of this problem, involving a numerical computation of the spectral function, 
is presented in Appendix A. This calculation confirms that the low-order perturbative result is qualitatively correct, 
even at large /o. Also in this Appendix we show how the Green's functions of the /o problem are obtained within the 
present generating function approach. Here, as a further illustration of this approach, we consider in detail a Fermi 
liquid with a /i interaction. The bosonized /i interaction contributes to the action the following term: 



SpL — 



/: 



2Vk 



S,T,q 



which can be separated into two parts, one longitudinal and the other transverse: 



S 



FL 



[kg ■ q] [kT • (-q)] , [kg x q] [ky x (-q) 



} J{T;-q) . 



(34) 



(35) 



Introducing two fictitious fields Ai and At which generate the interaction when integrated out of the partition function, 
the total action is given by: 



S[Au,a] 



E E 

S q,q-hg>0 



duj 
2^ 



{Lu-v*pci- ns) a*{S; q) a{S; q) + S'p^ 



(36) 



where 



S" — - 



1 



2/i 

1 



kp V 



duj 
2^ 
duj 
2^ 



d^q 



2 [Ai{q) Ai{-q) + At{q) Ai-q)] 



EE[^ M-q) 

S q 



ks X q 



J{S-q) Ati-q)] 



(37) 
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As usual, we contruct the generating functional and integrate out the boson fields, a and a*, to obtain the effective 
action for the An fields: 



where 



SefAAiA = IJ^J^ {[-!.+ ^,(q)] AM) M-q) + [-j-+xt{q)] A,{q) A{-q)} , (38) 



Xl{q) ^N*{Q)[~x^-l- + 0{x^ - l)\x\^^= - ie{l - x^)\x\- 



Xt 

The boson propagator is given by 

{a{S;q) a\S;q)) 



(q) = N*{0)[x^ - i - 9{x^ - l)\x\Vx^ - 1 - ie{l - x^)\xWl - x^] . (39) 



LO - v*pCi ■ fis + 177 sgn(i:j) 

. A ci-hs{q.--ns?Di{q) + {ci-KnsfDt{q) 



(27r)2 |q|2 [uj-v*pq,-ns + irisgTL{uj)Y 
where the propagators of the mediating fields are given by 



(40) 



These results will be used in the discussion of the physics of the half-filled Landau level given in the next section and 
in Appendix C. 

To complete this section we consider the effect of long-range, longitudinal interactions such as the Coulomb inter- 
action or the super long-range interaction of Bares and Wen.E3 Following the procedure given in Sec. we obtain 
the fermion Green's function by setting ifgQ = \/V{q) in Eq. ( |l8| ) and ehminating tho-jtransverse gauge fields. We 
briefly review what we find from this analysis and refer the reader to a previous papert3 for details. The correction 
to the boson propagator is given by 

' l + y(q)x*'(g) [^-i^ffis-q + «77Sgn(a;)]2 ^ > 

and the effect of screening is manifest. For example, in the case of the Coulomb interaction, V^(q) oc e^|q|^^-'^ and 
consequently, in dimensions greater than one, 5Gb is independent of the coupling constant in the limit of small 
q and small lo. We can calculate the correction to the fermion Green's function by expanding in powers of 5Gb- 
To first order, we obtain the well-known result that the imaginary part of the fermion self-euargy is proportional to 
[up' jep) ln|ijj| in D = 2, and {X/ephp) uj^ in D = 3, in agreement with Fermi liquid theorylHl. Below we will show 
that in the case of Coulomb interactions the weight of the quasiparticle pole, Zp, remains non-zero at the Fermi 
surface and the Landau Fermi liquid is stable as expected, ft might be thought that this is the case for all types of 
longitudinal interactions, as it appears to be ajjionsequence of screening. However, super long-range interactions can 
destroy the Fermi liquid state. Bares and WenEj considered a system of fermions in two dimensions interacting via a 
logarithmic potential, T^(q) = g/q^, and showed that within the random phase approximation (RPA) Zp — 0. Note 
that g is an energy scale and the plasmon gap is non-zero due to the super-long-range nature of the interaction and is 
given by ujp = vpy'mg/An. (By Kohn's theorem, the gap depends on the hare fermion mass and Fermi velocity. See 
Appendix C for details on how Kohn's theorem is automatically recovered within multidimensional bosonization.) 

One might be tempted to compute Zp by employing the Kramers-Kronig relation to derive the real part of the 
self-energy from the imaginary part estimated here. This procedure, however, is unreliable, as the Kramers-Kronig 
relations involve an integral over all frequencies, whereas the above calculation ignores high-energy processes such as 
the emission of plasmons by quasiparticles. Instead we compute the real-space, equal-time, two-point fermion Green's 
function directly. Because we wish to examine the effects of plasmons, we take the energy cutoff to be much larger 
than the plasmon energy, v*pX ^ Up. Obviously, then, bosonization is less accurate for systems with a large plasmon 
gap. Nevertheless, we expect it to be qualitatively correct, as discussed below. Setting t = 0, = x x fis = and 
introducing = x • 115, = q • n^, and (7_l = q x fig we find: 



8 



2:11 L 7_;,/2 27r 7-A/2 27r g2 

'■^ 1 1 



:} . (44) 



TT 1/5 — u|,m/(47rtj2) [uj — Vpq\\ + irj sgn{Lu)]' 

Here because plasmons are important only in the large-x limit, we have expanded the Lindhard function in powers 
of 1/x, X = w/(u|,|q|): 

=-^^ + 0(1/^4) +ir? (45) 

Kohn's theorem is respected by replacing the effective mass and velocity in the plasmon pole by bare quantities. 
Now since u) > Vp\q\, in the frequency integral we may make the approximation of replacing [uj — v^q^^ + irj sgn(w)]2 
by uj'^ + irj' and carry out the integrations. The integral over uj yields a constant due to the plasmon pole at non- 
zero frequency. Logarithmic dependence on X|| then is due to the factor of q'^ ^ q'^ + q'^^ in the denominator of the 
momentum integral. Upon exponentiating the boson Green's function we obtain the fermion quasiparticle propagator 
in configuration space: 

(A|a;|i|)-« 

Gf(5;x||) « e^'^^^ii ^ ' '"^ ; A|x||| » 1 
^11 

;A|x|||«l (46) 
•^11 

where 



The non-analytic dependence on a;|| demonstrates explicitly the destruction of the quasiparticle pole by the emission of 
plasmons. Power law behavior with exponent C, rather than the logarithmic behavior found in RPA, is a consequence 
of bosonization which treats the interaction non-perturbatively. In fact, Bares and WencJ used both a hydrodynamic 
calculation and a variational wavefunction to argue that the logarithm is promoted to a power law, with exponent 
C/V^, the difference by a factor of l/\/2 is due to spin degeneracy which is not taken into account here. This 
is further evidence that bosonization can accurately describe non-Fermi liquid systems. In the opposite limit of 
ujp\t\ ^ 1, however, an evaluation of the fermion propagator for a system of area yields the usual (a;|| — Vpt)~^ 
Fermi liquid form, multiplied by a prefactor of (AL)~^. As screening is effective at low frequencies, thermodynamic 
properties such as the specific heat which are defined in the equilibrium g-limit also retain the usual Fermi liquid form. 
Note that in either limit the propagator is odd under the combined CPT operation of (x, t) — s- — (x, t) and complex 
conjugation as demanded by Fermi statistics. The anomalous exponent C in the equal-time two-point function shows 
that the discontinuity in the quasiparticle occupancy at the Fermi surface has been replaced by a continuous, though 
non-analytic, crossover. The quasiparticle occupancy decreases from 1 to over a momentum scale of order A. Since 
the discontinuity has been smeared out, we might expect to encounter problems of consistency as the derivation of the 
current algebra Eq. (^) relies upon the existence of an abrupt change in the quasiparticle occupancy over a momentum 
scale much smaller than A. Nevertheless, as noted above, ±he form of the Green's function agrees with that found by 
an independent argument using a variational wavefunctionEj. The approximations implicit in bosonization (restriction 
to low energy processes, the existence of the discontinuity, and assumed stability of the fixed point) are robust in this 
case. 

It is instructive to repeat the calculation of the equal time Green's function for the case of a Coulomb interaction 
in both two and three dimensions. In two dimensions plasmons are found at a; oc ^/\q\. In this case the integral over 
to is no longer constant, as it was for the Bares and Wen interaction, but rather is proportional to \/|qJ and vanishes 
in the limit of small momentum. The momentum integral therefore no longer diverges logarithmically at large-xy; 

rather the integral is proportional to y^X/kp 0. While there is a non-zero gap for plasmon excitations in three 
dimensions, the integral over q again is suppressed, in this case because of the reduced phase space at small momentum 
and consequently SGb c>c X/kp 0. In either case the quasiparticle propagator has the standard Fermi liquid form. 
The weight of the pole undergoes no additional renormalization beyond the factor of Z]^ due to integrating out the 
high-energy degrees of freedom, Eq. (|^) . This is as expected, since the renormalization of the quasiparticle pole in 
Landau Fermi liquids is due to high-energy processes, not the low-energy processes found near the Fermi surface. 
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In this section we examined the effects of both short- and long-range longitudinal interactions on fermion liquids. 
We find that the Fermi liquid state is the only solution to the problem of a degenerate gas of fermions interacting 
via short-range or Coulomb two-body interaction in two and three spatial dimensions. Bosonization allows us to go 
beyond an assumed Fermi liquid form for the quasiparticle propagator. Indeed, had we considered fermions with spin, 
we would have found that the bosonizcd Hamiltonian separates into charge and spin partsllil, H — He -I- i?s, leading 
to the possibility that, as in one dimcnsioiiE3, the quasiparticle propagator would also exhibit spin-charge separation. 
Spin-charge separation in dimensions larger than one would, however, destroy the Fermi liquid, as the key element, 
the existence of a pole in the single-particle Green's function with non-zero spectral weight, would be replaced by 
a branch cut and the pole would be destroyed. This does not occur because the location of the pole of the boson 
propagator is unchanged from its free value in the A ^ limit. Consequently the sain and charge velocities are equal 
and both degrees of freedom propagate together in the usual quasiparticle formO. On the other hand, the Fermi 
liquid form is destroyed both in one dimension, where a Luttinger liquid occurs instead, and in the case of the super 
long-range interaction in two dimensions studied by Bares and Wen. The non-Fermi liquid fixed points of these two 
examples are accessible via bosonization. We therefore have reason to expect that bosonization can describe other 
non-Fermi liquid states of matter. The system we study in the next section is another such example. 



IV. TRANSVERSE GAUGE FIELDS: THE HALF-FILLED LANDAU LEVEL 



The two dimensional electron liquid in a magnetic field with even denominator filling fraction has been the subject 
of extensive studies. Experiments at even denominator fillings point_tp|-a, novel state which is compressible unlike 
the incompressible states found at odd denominator fractional fillings .EJEj Attempts have been made to understand 
the newjaauap-tum state by constructing composite particles which are made of magnetic flux tubes attached to each 
electrorO^tilEl. The composite particles see zero average magnetic field and in many respects behave like real particles. 
In this section we study the specific case of filling fraction v = 1/2. 



A. Half-Filled Landau Level 



The Hamiltonian for a two-dimensional electron gas of density in a perpendicular magnetic field B is given by: 

(48) 



H = I (fx cl{x)- i-^^^^ ^Ce(x) + Hcoulornb 



2m 

Here i? = V x A; at a value B = Airpfc/e half of the states in the first Landau level are filled. The HLR theoryS 
makes use of a local gauge transformation to describe this system as a collection of quasiparticles which obey Fermi 
statistics: 



cH^) = 4W exp ~i4> (fy arg(x - y) 4(y) Ce(y) 



(49) 



When (f) is an even integer the quasiparticles are fermions. Each quasiparticle is a composite object consisting of the 
physical electron together with a flux tube of integer quanta. The transformed Hamiltonian is given by 



H 



J-zV + (f)A-AO^ 
a X c'(x) ^ c(x) 



2m 



where 



A'(x) 



d-y- 



X - y) 



(x-y)2 



c^y) c(y) 



(50) 



(51) 



The average field strength of the flux tubes can now be adjusted to cancel out the external magnetic field. At 
half-filling setting (p = 2 



V X 



(-)A-A' 

c 



-'-iB 
c 



2-Kpfcq 



= 



(52) 



consequently, at the mean-field level, the quasiparticles behave as if they are free fermions in zero net field; the 
infinite Landau degeneracy is lifted and presumably the state is stable against two-body Coulomb interactions. Now 
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the important question is whether gauge fluctuations (or equivalently quasiparticle density fluctuations) modify or 
destroy this mean-field state. Unhke physical Maxwell electromagnetic gauge fields for which the gauge coupling is 
the fine structure constant a ~ 1/137, the coupling constant between the composite fermions and the transverse 
component of the statistical gauge field is of order unity. 

In the following we assume that despite strong gauge fluctuations the bosonic construction continues to apply. Later 
we will verify that this assumption is correct for the case of a Coulomb interaction. We determine the composite 
quasiparticle Green's function and show that the physics of the = 1/2 state is controlled by a marginal Fermi liquid 
fixed point. It should be noted that the gauge fields of the HLR theory arje-compact and consequently instantons iBiay 
play a role in the determination of the physical electron Green's functionO. The action of this system is given bylll: 

5-= I d^x dt \c*(i dt + Ao)c+ —c*(\7 + 1 (-)A-i A')^c+^(V x A')l 
J '- ' 2m ' 27r0 

- I J d^y dt [V. X A'(x)] F(x - y) [Vy x A'(y)] , (53) 

The last term in the action is the longitudinal two-body Coulomb interaction which now. is written in terms of the 
gauge field by making use of the constraint that a flux tube is attached to each electronEl, V x A' = — 27r0 c* c. It 
is important to notice that the longitudinal interaction plays a crucial role in the Chern-Simons gauge theory for it 
suppresses fermion density fluctuations. Since density fluctuations are equivalent to gauge fluctuations in a Chern- 
Simons theory, repulsive long-range longitudinal interactions moderate and control the gauge interactions. After 
integrating out the high energy quasiparticles following the procedure of Sec. ^ and adopting a change of notation 
A' — (e/c)A A, the HLR effective action reads: 

S= [ d^x dt [V/(i dt + Aq)^P + 7r^V'*(V - i A)V + ^(V x A)] 
J 2m* 211(1) 

^ ^ J <^^^ <^^V [^x X A(x)] - y) [Vy X A(y)] . (54) 



2 (27r0)2 

At filling fraction = 1/2 we set (p — 2 to cancel out the external field; now A should be interpreted as the fluctuation 
about this mean field value. Fermi liquid interactions should be included in the action also. If all coefficients except 
/i are set equal to zero, they are described by the term 



P S,T,q 



(55) 



which is made gauge covariant by the replacement^ k = — iV —iS/ — A. The bare gauge action, Eq.(|T^) can now 



be read off from Eq. (54) 



(56) 



here V{q) is the Fourier transform of V{x — y); for the Coulomb interaction, V^(q) = 27re2/|q|. For computational 
convenience, we work in the Coulomb gauge which eliminates the longitudinal component of the vector gauge field, 
only the transverse component A'^{q) = q x A((7)/|q| remains. 



B. Bosonization of Half-Filled Landau Level 



We now use the bosonization method to determine the composite quasiparticle two-point Green's function. As 
before, we bosonize the action Eq. ([s^), construct the generating functional and integrate out the boson 

fields to calculate the effective gauge action. The components of the inverse gauge propagator K^j^^ ~ [Z?^^]^^ (now 
the indices fi and v run over only the time and transverse directions, and T) are given by: 
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Kooiq) = X°{q) 

1 

Koriq) = --^|q| = Kroiq) 

ZTTip 

KMq) - (^^^(q) + X'^iq) ■ (57) 

Again the additive correction to the boson propagator is given by Eq. ( p6| ) where 6^(5*; q) is now the two-component 
vector (l,q x Vg/|q|). In the q-hmit, lo <^ the gauge propagator is dominated by the transverse components 

and 

here 7 = pfjixv *pN*{G) following the notation of HLR0. In this limit the gauge propagator has an imaginary part 
which is a result of quasiparticle damping. For the Coulomb interaction, 

B,M) e'iS;q) e'^iS;^q)^v* /'^''^'^' ^ (59) 

where x = 27re^/(47r^e0^), e is the dielectric constant, and the pole is located at w « i e'^q^ /kp- If, on the other h; 
V{q) is a short-range interaction, the pole is at uj i Vp |qp/fc^. The pole location agrees with the RPA resultl 
In the opposite a;-limit, u ^ ^^j^lql, we find 

D,A,) e^iS;q) eAS;^q) . i ^ ^^("^ ^^"^^^ ; (60) 

in this limit the poles which characterize collective excitationS, the cyclotron modes of the composite fermions, 
are found at = 2Tr pfcji/m*, the cyclotron frequency of a free particle of mass m* in an external magnetic field 
of i? = 2TTpfC(j)/e. As it stands this result violates Kohn's theoren£3. The correct result is obtained when the 
Fermi liquid interaction, Eq. (|5^), is included in the bosonized theory. The poles are then shifted to ^ ujc ~ 
2'Kpf{\/m* + /i/47r)(/) = 2'Kpf(j)lmh in accord with Kohn's theorem. The details of the proof of this result are given in 



Appendix C. The form of Eq.(60) is similar to that found for the super long range interaction of Bares and Wenc3 in 
the previous section. Since a projection onto the lowest Landau level is implicit in the HLR theory, however, we must 
assume v*p\ <C uJc so that inter-Landau level transitions arc not allowed. This restriction is an auxiliary condition 
which must be taken into account when defining the composite fermion Green's function. Non- Fermi liquid behavior, 
if it exists, must have its origin in thft behavior of the gauge propagator in the opposite g-limit. We may estimate 
the fermion quasiparticle self-energyO't3 from Eq. (^8|) by expanding the fermion Green's function in powers of D^y\ 
however it should be noted that there is no obvious small parameter to justify the expansion. Given this caveat, in 
the case of the Coulomb interaction at first order in the expansioijuwe obtain the marginal Fermi liquid (MFL) form 
for the self-energy in agreement with the RPA calculation of HLRfil. 

As the expansion in powers of the boson propagator is unreliable, it is more useful to determine the quasiparticle 
propagator directly in real space and time. We consider the fermion Green's function, Eq. (p^) and focus on the 
q-limit. We now set xi_ = and factorize the quasiparticle Green's function into the product of the free term 

G%{S-^,t)^ f'",, (61) 

and a factor which is the exponential of the additive correction to the boson propagator, SGB{S;:si,t), Eq. ([2^). 
Evaluation of the integral over g|| yields 

6GB{S;x,t) = 2i\x«\-^ sgn(x||)(x||/^,J-^) / - q± aq. 



Vp Jo Jj^ 91 - 



2 _ ,-|,.,|7 
X 



here 
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is a new momentum scale and 



C-V^ (64) 
is a dimensionless number. The final integration over frequency can be carried out using the integration formula 

/ e"^^ ln(/3 + x)dx^- [In /3 - e^'^Ei(-^/3)] (65) 
Jo 

where Ei(x) is the exponential integral function. The result is 

^^-(^^-'^)-- .||-.^.+lsgn(.||) {^ In[(.„ - .>.)^.^ + 4a^.^] 

_ g(.|i-«Jt)5 sg„(x||)+2^a^ _ (^^^ _ sgn(:E||) - 2ia5^ } (66) 

which, upon making use of the asymptotic expansion of the exponential integral function, reduces to 
5GB{S;^,t)KC\x\\5\ \n\x\\ ~ v*pt\5 ] \x\\ - v*pt\5 I 



-c 



x\ 



^11 ~ ^J"^ "f" ^^'^ sgn(a;||) 



"f + lii\xi\ - Vpt\S ; |a:|| - w^i|(5 > 1 . (67) 



Here 7 « 0.577216 is Euler's constant and a = 1/A. 

For simplicity, first consider the equal-time quasiparticle propagator given at long distances \ris • x| ^ 1/S by: 

GpiS; x) « 4^ (5 |ns • x|)-^ , (68) 
ns ■ x 

and at short distances jng • x| <C 1/(5 by the usual Fermi liquid form: 

iks-x 

Gf{S;x)^- . (69) 

ng • X 

At momentum scales less than d the quasiparticle pole is destroyed by transverse gauge interactions suggesting that 
the — 1/2 system is controlled by a novel non- Fermi liquid fixed point. Despite the appearance of the anomalous 
exponent C, the composite quasiparticle Green's function differs in an important way from the Green's function found 
for the super long-range longitudinal interaction in the previous section. The composite quasiparticle occupancy drops 
abruptly from one to zero at the Fermi surface over a momentum scale ^ oc which is much less than the cutoff A 
as A — > 0. As the basic consistency condition is satisfied in this limit, bosonization can provide a systematic approach 
to the study of the new fixed point. 

For typical GaAs/AlGaAs samples the carrier density pf is of order 10^^ /cm?, the electron band mass rrib — 0.068me, 
and the dielectric constant is given by e = 12.8. However, the effective mass of the composite fcrmion quasiparticle 
has to be estimated. Since electron states in the lowest Landau level are degenerate, the electron band mass does not 
play an important role. Instead, the energy scale is determined by the Coulomb energy at the magnetic length scale. 
From dimensional analysiaj the effective mass of the quasiparticle must be therefore: 



(70) 



e2 C 

where GJS-P' dimensionless constant to be determined cither numerically, or from experiment. Experiments at = 1/2 
suggestESE^ that m* « lOmf, and hence C « 0.10. For values in this range we find ^ = C0^/27r w 0.05. The small value 
of the exponent ^ is consister^with the abrupt change in the occupancy inferred from the experimental observation 
of sharp cyclotron resonancesEZIE2l. For the purpose of estimating 5, we may take the cutoff to be of order the Fermi 
momentum X ^ kp and find S — A^/(4C0^A;f) — 0{kp) as expected since kp is the only physical momentum 
scale in the problem. Note that the Fermi energy, ep — kp/{2m*), is still much smaller than the cyclotron energy, 
— kp(f)/{2mi,), and hence the requirement VpX << is respected even for A w kp. In the opposite limit |x| = 
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the quasiparticle Green's function is given by the free form GF{0,t) = l/{vp t). Consequently the density of states 
is the same as that of non-interacting fermions with Fermi velocity v*p. 

We examine briefly longitudinal interactions other than Coulomb. Because longitudinal interactions control fermion 
density fluctuations and hence, by the constraint, gauge fluctuations, we expect repulsive interactions which are longer- 
ranged than Coulomb to stabilize the Fcrmiliquid fixed point. For a general interaction ^(q) oc l/|q|^^^ this question 
has been addressed by Nayak and WilczckQ using a rcnormalization group expansion for small y. They concluded 
that for y positive (negative) the gauge interaction is relevant (irrelevant). The case y — Q is marginal and they 
found logarithmic corrections to Fermi liquid behavior. Here we set V{q) = 5|q|''~^ and repeat the calculation of the 
equal-time Green's function. We obtain: 



where x = 9 / i^iKf))^ . When y 0, the correction to the boson propagator scales as 5Gb oc {\x\\ \S')v/i^+v) where 
5' = ix/vp"/){vp/x)^^'^^^^y- Therefore, for y < 0, the gauge interaction is irrelevant and 5Gb — > as |a;||| oo, and 
the Landau fixed point is recovered. For y > 0, however, we expect pronounced deviations from the Landau fixed point 
when |(5Gb| S> 1, or equivalently when |a;||| ^ 1/(5'. The small momentum region is controlled by a non-Fermi liquid 
fixed point about which bosonization can tell us little as the deviations from Fermi liquid behavior are too ^large. 
For example, the discontinuity in the quasiparticle occupancy at the Fermi surface is completely eliminateclQ and 
inconsistencies develop within the formalism. The problem of a pv^e gauge interaction in two dimensions, considered 
by Lee and Nagaosa as a model for the cuprate superconductoral3, is the special case y — 1. Again this system is 
outside the scope of bosonization. 



C. A Marginal Fermi Liquid 

The real-space and time composite fermion Green's function given in the preceding section by Eqs. ( p^ , |6^ ) is the 
starting point for the more detailed discussion of single-particle properties which we give in this section. In particular 
we show that the Green's function is of a marginal Fermi liquid (MFL) type. In this section for simplicity we set 
Vp — 1 and denote X|| by a;. We first rewrite the Green's function in the two limits — i|(5 ^ 1 and \x — t\S <^ 1: 

GF{x,t)^ exp\C\x\d\n\x~t\S]; \x - t\6 <^ 1 

X — t + ia sgn(a::) 



X — t + ia sgn(x) x — t 



exp{-C ' ' }; |x-i|5>l. 



(72) 



Ideally we would Fourier transform the Green's function directly into momentum and frequency space to extract the 
self-energy. Unfortunately this task is quite difficult, so instead we begin with the observation that the pole of the free 



(5\: 



-t 



fermion Green's function at x = t has been destroyed by the transverse interactions. As a; — > t, Gpix^t) 
so the pole at x = t has been eliminated. To ascertain where the spectral weight is concentrated now, we expand 
the Green's function in the large-time limit \t\ ^ \x + (^x\n{\x — t\5)\ and show later that corrections arising from 
contributions at shorter times \t\ <C |a; -I- Ca;ln(|a; — t\6)\ contribute to the incoherent spectral weight. Since in the 
large-time limit \x/{x — t)\ <C 1 we may Taylor-series expand the exponential appearing in the second line of Eq. ^2 
Keeping the leading terms we obtain: 

G^(x,i)-( ^ ^ ^ 



x~t + ia sgn(x) ' i + S^^ \n{\x - t\S) + HOT 

= i f73) 

{x-t) + Cx\n{\x-t\5)+iasgn{x) + H0T ^ ' 

where HOT denotes higher-order terms which are negligible in the \t\ 3> |a;| limit. We may now Fourier transform 
this leading term to extract the self-energy. Upon making the approximation that ln(|a; — 1\5) w ln(|x|(S) near the pole 
in Eq. (u3), the integral over time can be done: 



GF{k,Lu) — J dxdt 



{x — t) + (x ln(|a; — t\5) + ia sgn(a;) 



14 



27ri / dx[ei-x)0{Lu) - e{x)e{-Lu)] e*('^-fe)^ exp[<tja;ln(|a;|(5)] 



2m ( dx[e{~x)e{u,) - e{x)e{-u;)] e^t-'^)- ['^^^^'^i'^l'^^]" . (74) 



ji=0 



Each term resulting from the Taylor series expansion of the exponential can be integrated separately. It is easy to 
show that: 

1"°° —1 

f dx e^'''-'^>[x\n{x5)f = jr^—^ {n! ln"[(77 - ifc)/(5] + 0(ln"-^[(r/ - zA:)/5])} (75) 



where > is a small term which ensures convergence. Up to subleading corrections, the Green's function is: 

LO — k 



n=0 



ln[i(cj — fc) sgn(ti>)/5] 



LO — k 
2tx 



{lu — k) — (^Lulii[i{LU — k) sgn(a;)/(5)] 
Hence the leading contribution to the self-energy of the quasiparticles is: 



(76) 



7 qr 

Y.{k,uj) = Cuj\ii{\Lo~k\/6) + —CoJSgn[{Lo-k)uj] . (77) 

At the pole, k ^ uj — ln(|a;|/(5), |a;| ^ |fc|, and to logarithmic accuracy, the self-energy is given by 

j:{k,uj) = CujH\Lj\/v*FS)--Cu; , (78) 

the marginal Fermi liquid form. In Eq. ( fz^ ) we have restored Vp. Note that Vp = kp/m* is the Fermi velocity 
obtained earlier by integrating out the high-energy degrees of freedom. In a marginal Fermi liquid such as this, there 
is no definite velocity; rather the velocity is a logarithmic function of the energy scale and is of order, but not equal 
to, Vp. At half- filling, (j> — 2 and ( = 2vp/Tre'^. Upon setting \ Ri kp for the purpose of estimating the size of the 
self-energy we find VpS ^ 0{£f), the Fermi energy; hence Eq. ( [ts] ) is identical to that derived by loffe et al. by a 1/N 
perturbation expansiontj. 

The evaluation of the self-energy might be questioned, in particular because the time integral in Eq. (Q) extends 
over the whole axis, whereas the approximate Green's function is accurate only at large time. Indeed, Gpixjt) does 
not have a pole at i = a; -I- Ca^ln(|a; — t\S) as can be seen from inspection of Eq. (p6|). To see why nevertheless the 
calculation is qualitatively correct, consider what happens if we restrict the time integral to large times by excluding 
intermediate times < t < 2tp for which the expansion of Eq. ( [73| ) is unreliable. Here tp = x + (^xhi(\x\S) and we 
can assume x > without loss of generality. The restricted integral over time may be evaluated giving: 



tp la J2tp ^ ^: 



{27ri6'(-cj) + Ei[(itp - a)uj] - Ei[-(«tp + a)uj]] 



= p('^-'tp)'^ [2Tiie{-uo) -Tri + 2iLutp + OiiLuhl)} (79) 



e ,^,,„>^, ^, ,, „ I ^„..^„„ , ^^^.^ Up/ 



where the second line follows from an expansion of the exponential integral about the origin Ei(z) =7-1- ln(2;) + z + 
z^/4 -I- 0{z^) with the assumption that tp ^ a. By neglecting the term proportional to 2iujtp and higher order terms 
which lead only to subleading corrections, the integral over x can be done. The resulting spectral weight is identical 
to that found from the imaginary part of Eq. ([76|). Apparently, coherent spectral weight in the {k,uj) plane is found 
in the long-time tail of the real-sgace and time MFL Green's function despite the absence of a pole in the {x, t) plane. 

In recent work Altshuler et alB have determined the quasiparticle propagator by a 1/iV perturbation theory in the 
fermion basis. In this work it is claimed that, because nonlinear terms in the fermion dispersion due to curvature of 
the Fermi surface are neglected in the bosonization scheme, the scheme cannot be used with confidence to study a 
transverse gauge theory. On the contrary, in the problem of the half-fUJed Landau level with a Coulomb interaction, 
the self-energy Eq. (Q) is identical to that found by Altshuler et ala Indeed, the only difference between the two 
calculations is the anomalous power law decay of the equal-time propagator discussed in the previous section which 
has its origin in the incoherent part of the quasiparticle spectrum and therefore should not be expected to appear in 
analyses such as that of Altshuler et al. which consider only the low frequency behavior of the Green's function. In 
the case of Coulomb interactions (y = 0) nonlinear terms in the fermion dispersion due to Fermi surface curvature do 
not affect the low-energy quasiparticle spectrum in any significant way. Why this should be the case is examined in 
the next subsection. 
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D. Consistency of Bosonization: Nonlinear Dispersion and Vertex Corrections 



An important feature of the bosonization procedure is the Unearization of the free fermion spectrum in the vicinity 
of the Fermi surface: eks+q ~ ^ks ~ • q. Nonhnear terms in the dispersion due-to curvature of the Fermi surface 
are neglected and within bosonization can only be accounted for perturbativelytj. This approximation might be 
questioned, particularly in the case of a mediating transverse gauge interaction, when the interaction vertex draws its 
most important contribution from q±, the momentum perpendicular to kg. To see this, consider the generic gauge 
propagator — i7|i^|/|q| — xIlT^^- The quadratic curvature term scales as q]^ oc and therefore it, rather than 
the pole familiar from the linearized approximation, v*pq\\^ might dominate the physics for y > 0. Thia^question is best 
addressed in the fermion basis by the consideration of the Ward Identities derived by Castellani, et ala. These authors 
found that the Green's function derived using the Ward Identities is identical to that found by bosonization provided 
that the relationship between density and current vertices, A5 = VgA^, which becomes exact in one dimension, holds 
approximately in dimensions greater than one. In the same context, we may ask also whether deviations from this 
relationship are so singular that they destroy the form of the propagator found both within bosonization and the Ward 
Identity approach. For the Chern-Simons gauge theory with y = 0, the physically relevant case of the long-range 
Coulomb interaction, we show that the bosonization scheme is internally self-consistent, and that the results derived 
earlier in this section continue to hold. 

We start by deriving the Ward Identities which relate the density and current vertices of the bosonized interacting 
fermion system in two dimensions; for simplicity, we give the proof for a fermion liquid interacting via a transverse 
gauge interaction only. The result is readily shown to be true in general. We introduce the time-ordered current- 
current correlation function between currents in different patches 

nS;Q(q) = -^(r J(5;q) J(g;-q)) (80) 

and the amputated vertex function 

r5;T(pi, P2; q) = -{T J(S; -q) ^(T; pi) V+(T; p2))amp (81) 
where q -I- Pi — P2 = 0. Then using the current algebra, we derive their equations of motion: 

- q ■ V* s)T s;TiPi, P2; q) = q - ns-^^ Vs-flXg) ^R;TiPi,P2;q) 

R 



and 



where 



+ 6s,t[g^\T;p2)~Gp\T;pi)\ (82) 



(i/ - q • w*g)Us,Qiq) = q • fis^ XI ^t-q{i) VQ;T{q) + Q ' i^s^^S,Q (83) 

T 



and Kj.rp[q) is the kernel of the bare gauge action. In addition, we make use of the Dyson equation which relates the 
reducible current vertices, Ts;T, and the irreducible current vertices, As-.t 

Ts;T{pi,P2;q) ^ As-T{pi,P2;q) + ^'ns:T{q) AR.T{pi,P2;q) Vq-Riq) . (85) 

R,Q 

In all the above equations, q which denotes (i/, q) is an index appearing in the bosonic quantities and pi denotes 
{LOi,Pi), an index for the fermion variables and for example, Gf{S;p) = G_F(a;,ks+p) with |ks| — kp ^ |p|. 
Combining Eqs. (b21 B3^, BSh gives the familiar Ward Identity relating the current and density vertices 



p A°{p;q)^ci-Asip;q)=Gp\S;p)-Gp\S;p~q) (86) 

where 
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(87) 



and 



As(p;g) = ^vJAt;s(p-9,p;9) 



(88) 



In addition, Eqs. (82, ^ ^5|) can be solved directly to determine Hs-.q^q): 



TT ( \ ^ ^'^S 



y J/ — q • vt 



17 q • ng q X vt, 

^S,Q + TT 



qx 



V ly-q-v* |q| /^^^^.(g) - p//m* - w^2^t(q) |q| 



(89) 



The diagonal elements of the correlation function are simply related to the boson Green's function Eq. (P^. These 
equations also provide an explicit calculation of the irreducible vertex of the bosonized interacting fermion system 



Gp\S-p2)-G-/{S-p,) 



from which it follows immediately that 

AT(p2;q) = ^v*g As;T{pi,P2;q) = v^A5^(p2;g) 



(90) 



(91) 



an exact result within bosonization but in general only approximately true. Introducing the quantity 

_ q- [As(p;g)-v* A%{p;q)] 



Ysip;q) 



mp;q) 



(92) 



and making use of the Dyson equation for the fermion self-energy and Eq. (|86|), Castellani, et alB found the fermion 
Green's function to be given by the solution of 



Gf{S]p) = GFn{S\p) + iGFo{S;p) 



dPq dv GpiS-jp — q) 



(27r)3 „-v*p qi\-Ys{p;q) 



D,,{q) e^{S;q) e''{S--q) 



(93) 



where G~^q{S\p) = o; — p • vj. When Ys{p; q) = 0, Eq. ( P3| ) can be solved and the result is exactly that derived by 
bosonization. 

Direct calculation of Ys{p;q) allows the investigation of the quality of approximation implicit in the bosonization 
approach. We evaluate Ys{p;q) in the fermion basis, to the lowest order, for a fermion liquid interacting via the 
Chern-Simons gauge field and a two-body longitudinal interaction of the form ^(q) = ff/lql^^^ and hence 



D^,{q) ef^iSiq) e-{S--q) 



2^txns| 



1 



■ kl I 



1+2/ 



(94) 



where 7 = pf /irvpN* {0} and x = 5/(27i'0)^. When the quadratic terms in the free particle dispersion are retained, 
the bare single-particle Green's function is given by Gp^{S;p) = — p • — p\/2m* . Correspondingly, the bare 
vertex function is given by 



A"5(p;'?) = i 
As(p;g) = vj 



2m* m* 



and therefore Ys{p] q) is nonzero at leading order 



yl{p;q) = - 



ql . p±q± 



2m* m* 



(95) 



(96) 



once curvature is accounted for. However, it is easily seen that these terms have no effect on the Green's function 
determined from Eq. (93). Integrating over q^^, Yg{p;q) exactly cancels the contributions due to curvature in the 
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Green's functions of the internal loop and the solution is as before. It remains to be shown whether or not vertex 
corrections due to the transverse gauge interaction modify this result. 

We follow Kim et alB and use a perturbative expansion, where N is the number of fermion flavors. The 
diagrams giving the correction to the vertex function at leading order in 1/7V are shown in Figs. 0, |^. At this order 
because the singular contributions arising from the two diagrams of Fig. |^ cancel each other by Furry's theorem, the 
only singular correction to Ys{p;q) is due to the diagram of Fig. |l| as can be verified by direct calculation. Furry's 
theorem states that for systems with chacge-conjugation symmetry, all diagrams containing a fermion loop with an 
odd number of vertices may be omittecO. Charge-conjugation (C) symmetry here is equivalent to particle-hole 
symmetry. The nonlinear pieces of the fermion dispersion, due to curvature of the Fermi surface, break this symmetry 
and the two diagrams of Fig. do not cancel completely. The remaining terms, however, are subleading and non- 
singular. From Eq. ( p3| , p6[ ), it is clear that the most important effect of finite Ys{p; q) on Gp comes from the pole at 
v — Vp g|| -|- q^/2m* — p±q±/m* = 0. At the pole, the density vertex is uncorrected and the only correction is to the 
transverse component of the current vertex: 



Asb; q)\ 



pole 



2m* 



— + n±5As±{p; g)|poic 
m* 



(97) 



where 



1 m*Vp 
q± "fTT^ 



A/2 



A/2 



-sgn(w) In ( 1 + 



dq' \q\ 

Hi 



sgn(..-.)ln(l + LII) 



q'\2+yX- 



(98) 



We see from Eq. ( [9^ ) that q±SAs± scales as (|'^|7/x)^^'-^^^''; substituted in Eq. ( |9^ ) this results in subleading 
corrections to the bare external fermion frequency, w, if y < and the fermion Green's function is unchanged. On 
the other hand, if y > 0, the corrections dominate and the results derived from bosonization are not reliable, as we 
had concluded for different reasons in Sec. 0. 

We now consider the marginal and physically interesting case of the longitudinal Coulomb interaction V{q) = 
27re^/|q|, y = 0, in detail. Carrying out the integration in Eq. (p8|), 



'5A5i(p;'Z)lpoic-^(4C)[c^lnf-z^) -(c.-i.)ln(-z''' ''^ 



v*p5 



v*pS 



subleading terms 



(99) 



substituting in Eq. (|93| ) and integrating over q|| we find 



Gf{S;p)^Gfo{S;p) 1 + 



{2n) 



dv I dq± 9[uuj) sgn{v) 



\UJ\ 



Yi{uj — v) 



(100) 



The q±_ integral over the gauge propagator is logarithmically singular at = and the z^-dependence in the denominator 
may be safely neglected on evaluating the final frequency integral. We then recover 



Gf{S;p)^Gfo{S-p) 1 + E(5;p) Gf(5;p) 



(101) 



where 



S-l \ ^1 f T-VpS + \UJ\\ 



(102) 



the marginal Fermi liquid form in Eq. (|7S 

In summary, for a fermion liquid interacting via the Chern-Simons gauge interaction and a long-range Coulomb 
interaction, the form of the propagator is not modified by quadratic terms in the quasiparticle dispersion due to curva- 
ture of the Fermi surface or by vertex corrections beyond those already incorporated in multidimensional bosonization. 
For this system bosonization provides a consistent calculational framework. 
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V. DENSITY RESPONSE FUNCTIONS 



In previous sections we focused attention exclusively on the single-particle properties of the interacting fermion 
system. In particular we found evidence of non-Fermi liquid behavior in the single-particle Green's function for the 
case of transverse gauge interactions. However, the single-particle Green's function does not have clear physical 
meaning as it is not gauge invariant. Indeed, as shown by Kim et al.a, different types of singularities might appear 
depending on the choice of gauge. These authors used a perturbative expansion to study such gauge invariant 
quantities as the density and current response at small momentum transfer. They were able to show that these 
properties of the interacting gauge system were indeed Fermi liquid like: singular corrections to the quasiparticle 
self-energy were canceled by singular vertex corrections. However, the possibility that singular behavior could arise at 
large manientum transfer was not ruled-put. In fact, other workers have found singular behavior at large momentum 
transferOE^. In particular, loffe et alSS have argued that the vertex correction is singular at momentum transfer 
2kp because the gauge fields mediate an attractive interaction between currents moving in the same direction. In 
this section, we calculate the gauge invariant response functions for both small momentum transfer Q ~ and for 
Q w 2fcj7. First we check that bosonization yields correct results for the free problem in both cases, and then we take 
into account the effects of interactions. 

We begin by writing the fermion density-density correlation function in bosonized form. The fermion density 
fluctuation is given by: 

<5p(x,t) = : V^t(x,t) V'(x,t) : 

= ^:Vt(S;x,i) ^(T;x,i) : (103) 

where the colons denote normal ordering which in this case is equivalent to the subtraction of the uniform fermion 
background density. Then the density-density correlation function can be written as: 



(5p(x, t) Sp{0, 0)) = ^Y.(J{S; X, t) J(T; 0, 0)) + {^f ^ ,^(ks-k.)x 



y2 "v-,","y/ J 

S,T S^T,U^V 



(104) 



For small momenta transfer, |Q| < A <C kp, and only the first t erm in Eg. (|104|) is relevant. The second term is 



important for momentum transfers near Q = 2kp. In subsections V A and VB[ we examine the first and the second 



terms respectively. We concentrate primarily on the case of two spatial dimensions. 

A. Density Response at Small-Q 

1. Free Fermions and Landau Fermi Liquids 

For free fermions, currents in different patches are uncorrelated. In two dimensions, and for |Q| <C kp-, the density 
response is given by the retarded correlation function 

HOo(Q,c.) = -z J dte^^--~-'^^Y.^J{S-^^) J(5;O,0)) 

= iV*(0)(^=f==-l), (105) 

where x = aj/(D|n|Q|). This is the well-known result for the density response of the free electron gas in the limit of 
low frequency and small momentum. 

Including a general longitudinal interaction V^(Q), we use the generating functional to calculate the correla- 
tion between currents in different patches. Then the retarded correlation function which combines contribu- 
tions from all possible two-point boson Green's functions (a(S';Q,w) a*(T;Q,w)), (a(5;Q,Ci;) a(T; — Q, — tj)), and 
{a*{S; Q, Lu) a*{T; Q, cu)) is given by 



19 



{J{S; Q, uj) JIT; -Q, -uj))n.t ^ ' ^ "^^^^ 



UJ ~ Vphs ■ Q + ir] 
A (lis • Q) {nx ■ Q) 



(27r)2 {lu - vpiis ■ Q + iri){uj - v*pnT ■ Q + ii]) 
V{Q) 



l-y(Q) n!]o(Q,w) 

Summing over S and T, we obtain 



identical to the RPA result for the density response which is exact in the limit Q,uj — > 0. 



(106) 



noo(Q,c.)^ wo! nolo V (107) 
1 - V{Q) n!jo(Q,w) 



3. Gauge Theory 



At two-loop order in a 1 /N pertubative expansion, Kim et a^Jl found that beyond RPA there were only subleading 
corrections to the irreducible polarizability. We showed in Sec. |IVD that implicit hi the bosonization scheme is a 
relation between the density and current vertices, Ag = vJA^. This relation suggestaj that within bosonization there 
are no corrections to the polarizability 11 beyond those of RPA. This supposition is correct as we will now demonstrate 
by explicit calculation for the case of the Chern-Simons gauge theory with a longitudinal interaction. 

The correlation function between currents in different patches can be found using the generating functional and is 
given by 



(J(5;Q,w) J(r;-Q,-a;))Ret 



in 
in 



ng ■ Q Ss^T 
oj - v*pns ■ Q + 
A 



[ns ■ Q) (fir • Q) 



(27r)2 [uo - v*p-n.s ■ Q + - VpUr -Q + iv) 



Q X v*g Q X 



IQI 



IQI 



Q 



IQI 



(108) 



where denotes the retarded part of the gauge propagator. Only the term proportional to Dqq{Q,uj) survives the 
sum over patches and we obtain for the density response: 



Uooici,Lo) = uUQ,^)[i + D^oiQ,^) nS]o(Q,^)] 



Qy{27r^y - nO,iQ,co) [pf/m* + Q^iQ) - nO,(Q, to)] 
where Il'j']^ (Q, cj) is the current response function of the free electron gas 

n°i (Q, c^) - ~v*p^N* (0) [x^-^-xVa 
In the same way, we may evaluate the current response function: 



1] 



where 



nn(Q,c^) = n°i(Q,c^)[l ~ Di^p{Q,iu) nUQ,uj)] 

nS]o(Q,w) 



n[!„(Q, Lu) [pf/m* + QWiQ) - n;,(Q, lu)] - QV{2TT~^y 
As stated above, we reproduce the RPA results, in agreement with the conclusion of Kim et al. 



(109) 



(110) 



(111) 



(112) 
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B. Density Response at Q ~ 2kF 



1. Free Fermions and Landau Fermi Liquids 



It is important to check whether bosonization captures the physics of the Kohn anomaly0. At first glance it is not 
obvious how bosonization, with the restriction |q| < A ^ fci?, can possibly describe processes with large momentum 
transfer, especially those near 2kp. The reason why it can is ma de cl ear by observing that scattering between 
different patches is allowed, and that in fact the second term in Eq. (104) contains these-processes explicitly. Recall 
that 2kp processes are described correctly by standard one-dimensional bosonizatioiH'CJ. Bosonization in one and 
higher spatial dimensions correctly describes processes with either small or large momentum, provided they are of low 
enerBMr-|That is why the specific heat, computed within the bosonization framework, is exact in the low-temperature 
limi#N. 

Since for free fermions the bosons in different patches are independent, the correlation function in the second term 
of Eq. (104) is zero except when V = S and U = T, and the Q m 2kp contribution to the density response function 
is given by: 



A 



D-l 



D 



ns • X — v*pt + irj sgn(i) 



A 



D-l 



i2n) 



E 



fix • X — Vpt H- ii] sgn(i) 



0{t) 



(113) 



Implicit in the last line is the restriction |x^A| < 1. As we are interested in the limit Q — ^ 2kF, only patches which 
are nearly opposite to each other on the Fermi surface participate in the scattering, hs ~ — fiT. Therefore we may 
carry out the sum over patches in each of the two bracketed terms in Eq. ( 113 ), keeping in mind that the spatial 
directions perpendicular to the Fermi surface normal are the same for both terms. To simplify the calculation, we 
consider only the static (w = 0) response and specialize to the case of a circular Fermi surface in two dimensions. 
Integrating the density response, Eq. (113), over positive time we obtain: 



n^,^(i?,c. = o) = 



'■^'^ ,^ ,^ exp[i/ci?i? (cos^i — cos6'2)][sgn(cos6'i) — sgn(cos6'2) 

d^id02 ^-—B 

f (costal — cos 02) n. 



(114) 



where 9i^2 are the angles between x and kg and respectively and R — |x|. Now we may focus on the case of 
Q w 2kp by limiting the range of integration ove r the angular variables 9i,2 to angles near 9i ~ and 02 ~ tt. In this 
limited region the cosines appearing in Eq. (114) may be expanded in Taylor-series as cos(6'i_2) = ±(1 — &\ 2/2) and 
the Gaussian integrals over the angles can be carried out. Finally, the Fourier transform into Q-space is carried out. 
If B is the angle between Q and x, the density response function is given by: 



no (Q,^ = o) 



2i kp-K^ 
wj^(2^)4 

1 J.1/2 



MR r,o'_^^kp^Qcose)R 
Jo 



8 v}.Tr3/^ 



in/4: 



dR 



i?2 

exp[i(2fcF - Q)R 



i?3/2 



(115) 



In the second line we have assumed Q >2kF and have also introduced an ultraviolet cutoff a to eliminate the artificial 
divergence at small- i?. The nonanalytic term we seek comes from the large- ii region defined by kpR ^ 1. Again in 
this limit a Taylor-series expansion of the cos 6 term is justified and the resulting Gaussian integral in the 9 variable 
has been carried out. For Q < 2kF, obviously Qcoa6 < 2kp and consequently the integral over 6 contains only 
terms which oscillate rapidly with increasing R and which therefore do not generate non-analytic behavior at large- i?. 
Finally, asymptotic analysis of the radial integral gives the non-analytic term in the density response function: 



^hkf 



(Q,^ = o) 



V2r 



Q-2kp 



47r V 2kp 
= C; Q <2kp . 



C; Q> 2kp 



(116) 



Here C is a constant which can be determined by a more careful evaluation of the integr als i n the small- i? region. 
Continuity across Q = 2kp ensures that the same constant appears in both limits of Eq. (116). The Kohn anomaly 
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of Eq. ( |116| ) agrees precisely with that found directly in the fermion basis.E3 Repeating the calculation in D = 3 we 
easily show that the static response is given by: 

n^.. (Q, - = 0) = ^ [ In + ^ _ 1 _ _] + C (117) 

which again agrees with the exact result found in the fermion basi^3 in the limit Q — > 2k p- The K ohn anomaly 
persists even if longitudinal interactions are included. In fact, this has already been anticipated in Eq. (|116|) and Eq. 



(117) where the effective mass m* rather than the bare mass m appears. 



2. Gauge Theory 

Now we consider the real-space density response function in two dimension at momentum transfer 2kp for the 
Chern-Simons theory. With the gauge interaction, currents moving in the same direction are attracted to each other 
so we might anticipate a singularity in 112^^^ . This singularity has a different physical origin than the Kohn anomaly 
of the previous section, which was not singular, just nonanalytic. The Kohn anomaly is due to the reduced availability 
of the phase space at low energies. In contrast, the singularity, if it exists, should be apparent even if we consider 
only scattering between two patches at exact opposite points of the Fermi surface. Therefore we need consider only 
terms with kg — — in the 2kp response function, Eq. (113). 



n2fcp(x,t) 



exp 



j- 47r 



[(^(5; X, t)(^{S- 0, 0)) + {<P{~S- X, t)ct){-S; 0, 0)) 



- (0(^;x,i)0(-5;O,O)) - (0(-^;x,i)0(5;O,O))] 

- ^ [terms with (x, t) (0, 0)] | e{t) . 



The in-patch boson correlation function is given by Eq. (E 



(118) 



{4>{S-, X, t)^(5; 0, 0) - 4>\S- 0, 0)) - ^ In 

47r 

+ -7-' 



X • ns - v*pt 



2^' 



2, |q X "si 



(119) 



The correlation function between bosons at opposite points of the Fermi surface contains the new physics. The 
correlation function {(f){S\ x, t)(f>{—S; 0, 0) — (j){S; 0, O)0(— 5; 0, 0)) can be found using the generating functional and is 
given by an integral over frequency and momenta of: 



g(g||) 

911 

n 9{q 



{e'(iii-ll—''> ~ l)(a(5;q) a{-S;-q)) + (e-*(9ll-ii-"*) - l){a\S;q) at(-5;-q))| 



+ (e-(«l,-,i-*) _ 1) D^,{q) - v*p^{D^^{q) + DUq)} + 



Vp qj 



(120) 



As in the case of th e two-point function, the most singular contribution comes from the g-limit. Inserting Eqs. (HE, 
|120|) into Eq. (118) and computing the integrals in g- limit we obtain 



U2k^ (x, i) = XI T ' TTT" } ^TTl ' — , , ■ i ^ ^^^P ' , 0] 

^ |_x • ns - v*pt + 17] sgn(x • 115)] |_x • ns + v*pt + i-q sgn(x • ns)J 



cx 



(|x| - v*pt + ili) (|x| + v*pt + iyy) 



exp [£(|x|,i)] , 



(121) 
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where the prime on the above sum over patches indicates that the sum is only over patches for which |x x figjA < 1 
and where 

£{x, t) « -C| ^ In I {x - v*pt)5\ + — In |(x + v*pt)6\ 

- In^ I (x - v*pt)5\ ~ \t? I (x + v*pt)5\ | , (122) 

for — v*pt\5 ^ 1 and jxy + v*pt\5 3> 1. The log-squared terms appearing in £ have their origin in the correlations 
induced by the gauge interaction between the bosons in opposite hemispheres. Evidently a log-squared singularity 
at momentum transfer 2kp exists in the real-space response function. Although it is technically difficult to Fourier 
transform it into frequency space, we can estimate the leading singular factor of 112^^ (w) with the use of a scaling 
argument and find: 

n2fc^(c^) oc exp(2Cln2 \uj/{v*p5)\) (123) 
where again v*p5 — 0{ep). This singularity is similar to, though stronger than, that found by Altshuler et a/M 

VI. TRANSVERSE GAUGE FIELDS: FERMIONS IN THREE DIMENSIONS 

In this section, we discuss the most familiar physical example of fermions interacting via gauge interactions: electrons 
in three spatial dimensions with the full Maxwell electromagnetic field. The longitudinal Coulomb interaction is 
screened, and the transverse gauge interaction is suppressed by a factor of fine structure constant multiplied by the 
ratio of the Fermi velocity to the speed of light. For these reasons we expect the behavior to be controlled by the 
Landau Fermi liquid fixed point except at extremely low temperatures. Nevertheless it is an interesting exercise to 
study the crossover to the non- Fermi liquid fixed point. l— , 

An expansion in the fine structure constant leads to a modified quasiparticle spectrumE3 near the Fermi surface: 

e{v)^VFM-pF)-ae^{vF/cfM-PF)\^\\p\-PF\ , (124) 

here a is -a. positive constant. It has been conjectured that the logarithmic term exponentiates to give a power law 
spectrumo, 

e(p) = vf{\p\ -Pf)\\p\ ^PFr"""'^' ■ (125) 

Here we investigate the fixed point using the formalism developed in previous sections. However, we do not find a 
power law, rather, a logarithmic spectrum of marginal Fermi liquid type characterizes the fixed point. 
We consider the following bare gauge action in the Coulomb gauge: 



1 



87r(e/c) 



1 ,dA 



Sg[a^ = iizrrr^2 / ^ (^)' - >< ' (126) 



where the longitudinal interaction is neglected since it is screenedQ. Once the bare gauge action, Eq. (p^ , is specified 
we can follow the prescription given in Sec. || to obtain: 



|q X nsp A'n:v*pe^ 

|qp 47re^ x'^((j') + cj^ — c^lqp 



where the transverse susceptibility il) was defined in Eq. (EG) . Th is re sult is in agreement with the RPA calculation 



of the gauge propagatoi£3. In cj-hmit, uj > v^|q|, the poles of Eq. ( |127D are at = Airpfe^ /m + c^|qp, the plasma 
frequency. In the opposite q-limit, uj < Vp\q\, we obtain: 



D,.(g) e^S) e''(S) = ^ ''^iJ'''. , , (128) 



where now 7 — 7rvpN*{0)/4 and x = c^/47re^. Inserting this expression into the boson propagator Eq.(p2[), we obtain 
the imaginary part of the boson self-energy due to quasiparticle damping. As in the case of the half-filled Landau 
level, the g- limit is the more singular of the two limits, and is the origin of the non- Fermi liquid fixed point. From 
Eqs.(p2l m, |12^), we can now compute the real space fermion Green's function. Factorizing the quasiparticle Green's 
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function into a product of the free term (S; x, = (x • fis — Vp t) ^ and the correction term exp[^GB(S; x, t)], we 



have 



(5GB(S;x,t) = i 



(2^)2 
(2^)2 



dq\\ 
2tt 



duj 
2^ 



?(q-x— 



11 



3c2 



27r 
2tt 



J^.(q,c.) e^(S;g) e'-(S;-g) 
[cj ~ q • lis + i77 sgn(w)]2 

iui sgn(a:ii) {x«/vp-t) 



e 

ii^ sgn(a;||) (xy/Dj,-* 



In 



x\<i±r - 



-lUJ 



(129) 



in the second line we use the fact that |q_L | 3> \q\\ \ in the region of integration where the integrand is the most singular, 
and perform a contour integration around the pole located at g|| — u!/v*p. The frequency integration in the last line of 
Eq.( |l29[ ) is identical in form to Eq. (|62|). Introducing a new momentum scale 5 = \^ / Ae'^v*pk'^p and dimensionless 
constant C = 2e^u|,/3c^, we obtain: 



(5GB(S;x,i) « -C|fe||| hi\5{x\\ - v*pt)\ ] \x\\ - v*pt\5 <Si.l 



< 



x\\ — v*pt + 2ia sgn(a;|| ) 



7 + In |(x|| — i;^i)(5 sgn(x||)| ; |a;|| — u^t|(5 > 1 



(130) 



As in the case of the half-filled Landau level, the Green's function has the marginal Fermi liquid form with self-enerejj 
given by Eq. (|78|). A similar result was found in a fermion renormalization group calculation by Chakravarty et. als^ 
However, because C, ~ 10~^ for a typical metal, the crossover to MFL behavior occurs at an exceedingly low energy 
scale Lo set by 



C \og{E/oj) = 1 



(131) 



where E = (c^kp)/ {Ae'^). Including nonlinear terms due to curvature of the Fermi surface in the fermion disper sion 



and vertex corrections due to the gauge interactions, we may apply the arguments presented in Sec. IV D. It 



is straightforward task to show that bosonization is consistent and provides an adequate description of the three 
dimensional gauge- fermion liquid. 



VII. CONCLUSIONS 



In multidimensional bosonization, Feynman diagrams are replaced by geometrical and algebraic ways of under- 
standing strongly correlated fermions. The generating functional approach is a convenient way to compute correlation 
functions of bosons and hence, fermions, interacting via short- and long-range longitudinal, and transverse gauge 
interactions. In one dimension well-known results for the Luttinger liquid are recovered. Short-range and Coulomb 
interactions in two or higher spatial dimensions do not destroy the Landau Fermi liquid fixed point; spin-charge 
separation does not occur. The Kohn anomaly in the density response function near momentum 2k p has the standard 
nonanalytic form, showing that all low energy processes, even those with large-momentum, are correctly described by 
multidimensional bosonization. Likewise, the specific heat has the correct linear dependence on temperature with sub- 
leading corrections which are easily computed.Ej Collective modes in the chacge and spin sectors follow directly from 
the equations of motion for the abelian U{1) and non-abelian SU{2) currentsO The stability of the Landau fixed point 
against superconducting BCS, CDW, and SDW instabilities also can be ascertained using the renormalization-group 
in the boson basisO Novel fixed points in higher dimension arise in several cases. These examples demonstrate that 
bosonization is more encompassing than Landau Fermi liquid theory as non-Fermi liquid fixed points are accessible. 
For a super-long range longitudinal interaction in two dimensions, a system first studied by Bares and Wen,L3 we find 
an anomalous exponent in the equal-time Green's function. For transverse gauge interactions we find marginal Fermi 
liquid (MFL) fixed points in two instances: for the HLR theory of the half-filled Landau level with the Coulomb 
interaction and for ordinary Maxwell gauge interactions in three dimensions. 

We studied the gauge-invariant density response function at the half- filled Landau level fixed point. At low momen- 
tum and frequencies we obtain the same result as the RPA calculation, and the response function is indistinguishable 
from that of a Landau Fermi liquid. Deviations from Fermi liquid behavior appiear near momentum transfers of 2kp. 
Since odd-denominator fractional fillings are equivalent, via Jain's construction,c3 to the half-filled Landau level in an 
additional external magnetic fiekl AB, it would be interesting to examine questions of the divergent renormalization 
of the composite fermion massEsE^ and the associated compressibilit}C3 as a function of AB within the bosonization 
framework. 
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The intimate connection between multidimensional bosonization and the Ward Identity method of Castellani, Di 
Castro, and Metzner enables us to go outside the framework of bosonization to check its consistency and accuracy. 
Neither non-linear terms in the fermion dispersion due to Fermi surface curvature, nor vertex corrections due to gauge 
fields, change our results in any qualitative way. All of the fixed points we found possess much higher symmetry than 
the bare interacting problem. This infinite C/(l) symmetry reflects vanishingly small large-angle scattering at such 
zero temperature quantum critical fixed points. Multidimensional bosonization is a good description of Fermi and 
nearly- Fermi liquids. 
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APPENDIX A: NUMERICAL COMPUTATION OF THE SPECTRAL FUNCTION 



Small-angle scattering between the fermion quasiparticles generates a subleading additive correction to the free 
boson propagator. We use the bosonization formula Eq. ( p^ ) to compute the fermion spectral function. Since 
bosonization is carried out in (x, t) space we must perform three operations. First we inverse Fourier transform 
the boson propagator into real space and time. Next the exponential of the resulting expression yields the fermion 
propagator. Finally a second Fourier transform of the fermion propagator back into momentum space allows us to 
extract the spectral function. 

It is a difficult technical problem to carry out these three steps analytically. In this Appendix we resort to the 
numerical method of fast Fourier transforms (FFTs) to compute the full spectral function. The problem is still 
challenginffi-since we seek to extract sub-leading corrections to scaling. For example, the imaginary part of the fermion 
self-energyEJ is of order {uj'^ /ep) ln(|a;|) and therefore much smaller than the quasiparticle energy at small frequencies. 
In fact, for spectral broadening to show up in the numerical calculation, it is necessary to set A — 0{kF)- As we shall 
see, setting A = 0{kF) forces us to treat velocity renormalization carefully. 

As a illustrative example of the numerical method, we consider the problem of a two-dimensional spinless Fermi 
liquid with only a single non-zero Landau parameter, /q. In this case we keep only the temporal component of the 
gauge fields appearing in Eq. ([l^) and set K^f^{q) = 1/fo- Then the gauge propagator appearing in Eq. ( [2^ ) is given 
by: 



A 

l + /oX°(9) 



D^.(q, Lo) e^^iS; q) e^{S- -q) = ^ , . (Al) 



To continue our analysis, we now set x± = in Eq. ( p4| ) as the spectral function has only weak and uninteresting 
dependence on transverse momentum. With this simplification, the integral over q± in Eq. (|2^) may bej^ierformed 
analytically, yielding a rather complicated expression which, for the sake of brevity, we do not present hereo. Instead 
we denote the integral by: 

Diqn,^)^ ^ (A2) 

which, multiplied by the double pole [uj — v'p q\\ +iri sgn(ijj)]~^, is the quantity to be Fourier transformed into {x\\,t) 
space. Here we distinguish between v'p, which is the velocity of the non-interacting bosons, and Vp, which is the final 
velocity of the bosons (and fermions) after the effect of the current-current interaction /o on the velocity has been 
accounted for. 

For Landau Fermi liquids, the fermion velocity renornializes from its free value when high-energy processes, those 
with energy greater than the scale set by the cutoff A, are integrated out. As long as A <C kp, further renormalization 
of the velocity due to the remaining low-energy processes, the ones to be bosonized, is insignificantO. Velocity 
renormalization cannot, however, be neglected when A is comparable to the Fermi momentum. To incorporate this 
renormalization directly into the Green's function, from -D(9|| , t^) we subtract its value at q\\ = lo = 0, D(0, 0). Setting 
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Fq = /oiV'(O), where iV'(O) = kp/^nvp, we find i?(0,0) ^ j^^. The final renormahzcd Fermi velocity is then 
given by: 

"F = (1 + ^T^)«F- (A3) 
2,T:kp 1 + _f 

The diff'erence Z)((7||,a;) — -D(0,0), which is small for small frequencies and momenta, now appears in the numerator 
of the remaining integral over u) and g|| in Eq. (p^: 



These two integrals may now be performed numerically by a two-dimensional (inverse) FFT. In practice it is simplest 
to study either the advanced or retarded Green's function rather than the full time-ordered Green's function; hence 
we restrict the FFT to the half plane t > Q. Upon subtracting the x\\ = t = value, exponentiating the difference, 
and multiplying by the free Green's function, the full fermion Green's function is obtained in real space and time. 
This result may then be conveiied back to {k,uj) space by the second FFT. The spectral function is then extracted 
by taking the imaginary parto. Displayed in Fig. || are cross-sectional plots of the spectral function at different 
frequencies computed on a 2048 x 2048 lattice. The peak has been centered in each plot. The width of the peak at 
zero frequency is due solely to finite-size effects. As the frequency increases from zero, the width grows rapidly as 
expected. The nonperturbative results demonstrate that low-order perturbative expansions in the Landau parameters 
are qualitatively correct as there is no breakdown in Landau Fermi liquid behavior. 

APPENDIX B: GENERAL LANDAU PARAMETERS IN TWO AND THREE DIMENSIONS 

It is straightforward to include higher order Landau parameters in the computation of the 2N-point boson correlation 
functions. Extension to three spatial dimensions is straightforward. Finally, Landau parameters in the spin sector may 
be incorporated in a perturbative way with the use of Abelian bosonization. In two dimensions the boson propagator 
for a single /„ interaction in the charge sector is given by: 

{a{S-q) a^{S-q)) 



where 



LO — v*p(\ ■ hg + irj sgn(aj) 
. A . cos^jn (3{S;q)]D[''\q)+sm-'[n /3(5; q)]!),^"^ (g) 

\27r)2'l'"^ [c. - «J,q • fis + sgnH]2 ' ^ ' 

D["\q) ^ , , , (B2) 

-!//« + Xi (<?) 

D['"\q) ^ , , , (B3) 

-l/fn + x't'\<l) 



and 



dd cos 9 cos^(n 



2n X — cos 9 + ir] sgn(x) 

,2^ 



X^\q)^N*iO)l'^ COS. sin (n.) 

' 271 X — COS 9 + IT] sgn(x) 



The function /?(S; q) in Eq. (Bl) denotes the angle between S and q. In three dimension, the inclusion of ^th Landau 



interaction can be performed using the addition theorem. The boson propagator is then: 

(a(S;g) at(S;g)> = , — 

uj - Vpq -ns + irj sgn(Lj) 

,.2A2 q.fts ^ Yr"'(^^(S;q)) Yr(n(S;q)) 
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where 



X'^^q) = / d^n ^ Y-™(n) Y^m (B6) 

' ^ ' An J a; - cos 6* + ?ry sgn(x) i ^ ^ « ^ ^ ^ ^ 

Here, r2(S;q) denotes the angles between S and q in spherical coordinates. 

In general, when more than one Landau parameter is involved, we introduce the appropriate number of mediating 
fields for each Landau parameter. For example, in two dimensions we introduce two mediating fields for a single 
Landau interaction /„ as above. In three dimension, on the other hand, we decompose the Landau interaction as 

(2/ + 1) fi Pi{coae)^An fi J] (-1)™ Yr(l1(S;q)) Yr™(n(S;q)) (B7) 

m— — l 

and introduce a mediating field A'^{q) for each spherical harmonic, amounting to 2Z + 1 fields for each Landau 
interaction. The propagators of the mediating fields are then matrices. Upon the inclusion of gauge interactions along 
with the Fermi Liquid interactions, we first construct a gauge- covariant form for the Fermi liquid interaction and 
then treat the gauge fields and the mediating fields on the same footing. An illustrative calculation is presented in 
Appendix C for the case of the Chern- Simons gauge interaction plus the Fermi Liquid parameter /i. 



APPENDIX C: HOW KOHN'S THEOREM IS SATISFIED 



In this appendix, we prove the result quoted in Sec. IV B that the collective excitations of the interacting quasi- 



particles of the half-filled Landau level occur at a frequency lo = 2T:pf(j)/m, the cyclotron frequency of a particle with 
bare mass m in a magnetic field B — inpfccj)/ e. To prove this, we add to the action of the Chern-Simons gauge theory 
a Fermi liquid interaction in which all the coefficients except /i have been set equal to zero. The contribution to the 
action Eq. ( ^5|) is given in real space and time by 

Sp^^J^ [ (fx dt tl^^Vtlj ■ ij^V^P (CI) 

which is made gauge covariant by replacing V by V — lA. On Fourier transforming it is given in terms of the currents 
as 



+ ^ k5-A(-g) + J(5;-g) k5-A(g)] Y,p) Hi) ■ M~q)] (C2) 



As in Sec. Ill the term quadratic in J is decoupled by introducing additional transverse At and longitudinal Ai 
fictitious fields. 

Including the Fermi liquid interaction the complete action of the theory is given by 



5[A^a,C,r] =E E {{u^v*pci-ns) a*{S;q) a{S; q) + ^S; q) a*{S;q) + C{S;q) a{S;q)} 

S q 

, n5 -q ns X q i 

+ ±y f^\[PL + ll^ + SliCMl Ar{q) Ari^q) + Aoiq) Ari-q) 

^2V^J 27r Um* kj, (2^0)2 J ' ^' 27r0 ^ ' ^ ' 

-^Ai {q) Ai {-q)-^ At {q) At (-g) } , (C3) 

where fj, runs over 0, T, I and t. and we couple the boson fields to ^, ^* to construct a generating functional as in 
Eq.(n6h. After integrating out the boson fields a, a*, we obtain an effective action: 



27 



Xtiq)] A{q) Ai-q) 



here 



+ ^Ao{q) Ar{-q) + {-^ 

ZTTCp Jl 



Xl{q)] Mq) M-q) 



+ xioiq) [Mq) Aa{-q) - Ao{q) Ai{^q)] + 2[v*p 
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(C4) 



(C5) 



and the susceptibilities Xa{Q) cire defined by: 

rif) 

X'iq)^-N*i0) ■ 



2tt X — cos 9 + irj sgn(a;) 



= N*{0)[l-e{x^ -1) 



Xt{q) = N*iO) 



2 / 
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27r a; — cos 6 + irj sgn(a:) 
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27r a; — cos 6 + irj sgii(a;) 
N*iO)[-x'-l + 0{x'-l)\x\^^ 
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d6i 



2-K X — cosO + 177 sgn(x) 



iV*(0)[-a; + 6i(a;2 - l)|x|- 



and X — LD / {v*p\(\\) . The boson propagator is obtained in the usual way 



ie{i - x'^)\x\ 
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(C6) 



(C7) 



and the collective modes ar e given by the poles of the gauge propagator. The inverse gauge propagator can be read 
off from the action Eq. (C4) and in the limit of interest ^ z;|,|q| is given by 



[K{q)U 
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fl 



•m*_hPf\ 



(C8) 



The poles of Dct/3 are found by setting det|ifQ^| = which gives cj^ = [2Tr(f)pf{l/m* =t-,/i/47r)]2 — {2t: p f (j) / . The 
second equality follows from the application of Galilean invariance to a Fermi liquidEJ. The collective mode appears 
at the cyclotron frequency determined by the bare electron band mass in agreement with Kohn's theorem. 
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A copy of the "C"-code executing these operations is available from the authors. 
FIG. 1. The first of two 1/N vertex correction, A''^^(p; q). The symbol X indicates that external lines are to be amputated. 

FIG. 2. The two diagrams which give the second 1/N vertex corrections, A^^^(p;q). The leading divergences in these two 
diagrams cancel via Furry's theorem (see text) as only the nonlinear terms in the fermion dispersion break charge-conjugation 
symmetry. 

FIG. 3. Spectral function lmG{k,u!) in two dimensions for the case _Fo = 2. Slices at three different frequencies w = 0,^^/8, 
and A;f/4 are shown respectively in (a), (b) and (c). Here A = kp, Vp = 1, and a lattice of 2048 x 2048 points is used to 
compute the Fourier transforms. The peak has been centered in each plot around the point k = uj. 
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